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0^ I Abstract. In these lectures, we discuss different types of renormalization problems 

0^ ■ in QCD and their non-perturbative solution in the framework of the lattice formula- 

tion. In particular the recursive finite size methods to compute the scale-dependence of 
i renormalized quantities is explained. An important ingredient in the practical applica- 

Q ■ tions is the Schrodinger functional. It is introduced and its renormalization properties 

' are discussed. 

Concerning applications, the computation of the running coupling and the running 
quark mass are covered in detail and it is shown how the yl-parameter and renormal- 
ization group invariant quark mass can be obtained. Further topics are the renormal- 
' , ization of isovector currents and non-perturbative Symanzik improvement. 
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1 Introduction 

The topic of these lectures is the computation of properties of particles that are 
bound by the strong interaction or more generally interact strongly. The strong 

1 Lectures given at the 36. Internationale Universitatswochen fur Kern- und Teilchen- 
physik (Schladming 1997) : Computing Particles 
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interactions are theoretically described by Quantum Chromo Dynamics (QCD), 
a local quantum field theory. 

Starting from the Lagrangian of a field theory predictions for cross sections 
and other observables are usually made by applying renormalized perturbation 
theory the expansion in terms of the (running) couplings of the theory. While this 
expansion is well controlled as far as electroweak interactions are concerned, its 
application in QCD is limited to high energy processes where the QCD coupling, 
a, is sufficiently small. In general - and in particular for the calculation of bound 
state properties - a non-perturbative solution of the theory is required. 

The only method that is known to address this problem is the numerical 
simulation of the Euclidean path integral of QCD on a space-time lattice. By 
"solution of the theory" we here mean that one poses a well defined question like 
"what is the value of the 7r decay constant", and obtains the answer (within a 
certain precision) through a series of Monte Carlo (MC) simulations. This then 
allows to test the agreement of theory and experiment on the one hand and 
helps in the determination of Standard Model parameters from experiments on 
the other hand. 

Quantum field theories are defined by first formulating them in a rcgulariza- 
tion with an ultraviolet cutoff yl cut and then considering the limit A cut — > oo. In 
the lattice formulation (Wilson 1974), the cutoff is given by the inverse of the 
lattice spacing a; we have to consider the continuum limit a — > 0. At a finite 
value of a, the theory is defined in terms of the bare coupling constant, bare 
masses and bare fields. Before making predictions for experimental observables 
(or more generally for observables that have a well defined continuum limit) the 
coupling, masses and fields have to be renormalized. This is the subject of my 
lectures. 

Rcnormalization is an ultraviolet phenomenon with relevant momentum scales 
of order o _1 . Since a becomes weak in the ultraviolet, one expects to be able 
to perform renormalizations perturbatively, i.e. computed in a power series in 
a as one approaches the continuum limit a — > 0. 2 However, one has to take 
care about the following point. In order to keep the numerical effort of a simu- 
lation tractable, the number of degrees of freedom in the simulation may not be 
excessively large. This means that the lattice spacing a can not be taken very 
much smaller than the relevant physical length scales of the observable that is 
considered. Consequently the momentum scale a -1 that is relevant for the renor- 
malization is not always large enough to justify the truncation of the perturbative 
series. In order to obtain a truly non-perturbative answer, the renormalizations 
have to be performed non-perturbatively. 

Depending on the observable, the necessary renormalizations arc of differ- 
ent nature. I will use this introduction to point out the different types and in 
particular explain the problem that occurs in a non-perturbative treatment of 
renormalization. 

2 For simplicity we ignore here the cases of mixing of a given operator with operators 
of lower dimension where this statement does not hold. 
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1.1 Basic renormalization: hadron spectrum 

At this school, the calculation of the hadron spectrum is covered in detail in 
the lectures of Don Wcingarten (Weingarten 1997). I mention it anyway be- 
cause I want to make the conceptual point that it can be considered as a non- 
perturbative renormalization. I refer the reader to Wcingarten's lectures both 
for details in such calculations and for an introduction to the basics of lattice 
QCD. 

The calculation starts by choosing certain values for the bare coupling, g , 
and the bare masses of the quarks in units of the lattice spacing, am„. The flavor 
index / assumes values / = u, d, s, c, b for the up, down, charm and bottom 
quarks that are sufficient to describe hadrons of up to a few GeV masses. We 
neglect isospin breaking and take the light quarks to be degenerate, mJJ = = 
m l . 

Next, from MC simulations of suitable correlation functions, one computes 
masses of five different hadrons H , e.g. H — p, 7r, K, D, B for the proton, the pion 
and the K-,D- and B-mesons, 

amjj =amif(jo,amo,aWg,omQ,am[5) . (1) 

The theory is renormalized by first setting m p = ?rip Xp , where m° xp is the exper- 
imental value of the proton mass. This determines the lattice spacing via 

a = (am p )/m p xp . (2) 

Next one must choose the parameters am^ such that (1) is indeed satisfied with 
the experimental values of the meson masses. Equivalently, one may say that at 
a given value of go one fixes the bare quark masses from the condition 

(am H )/(am p ) = m^/m^ , H = n, K, D, B . (3) 

and the bare coupling g then determines the value of the lattice spacing through 
(2). 

After this renormalization, namely the elimination of the bare parameters in 
favor of physical observables, the theory is completely defined and predictions 
may be made. E.g. the mass of the Z\-resonance can be determined, 

m A = a-^am^ [1 + 0(a)] . (4) 

For the rest of this section, I assume that the bare parameters have been elimi- 
nated and consider the additional renormalizations of more complicated observ- 
ables. 
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Note. Renormalization as described here is done without any reference to pertur- 
bation theory. One could in principle use the perturbative formula for (a A) (g ) 
for the renormalization of the bare coupling, where A denotes the /1-parameter 
of the theory. Proceeding in this way, one obtains a further prediction namely 
m p /A but at the price of introducing 0(<?q) errors in the prediction of the observ- 
ables. As mentioned before, such errors decrease very slowly as one performs the 
continuum limit. A better method to compute the yl-parameter will be discussed 
later. 

1.2 Finite renormalization: (semi-)leptonic decays 

Semileptonic weak decays of hadrons such as K — > ir e v are mediated by elec- 
troweak vector bosons. These couple to quarks through linear combinations of 
vector and axial vector flavor currents. Treating the electroweak interactions at 
lowest order, the decay rates are given in terms of QCD matrix elements of these 
currents. For simplicity we consider only two flavors; an application is then the 
computation of the pion decay constant describing the leptonic decay tt — > e v? 
The currents are 

A^x) =^{x) W \t^{x) , 

W=^(*)7„§t»#c) , (5) 

where r a denote the Pauli matrices which act on the flavor indices of the quark 
fields. A priori the bare currents (5) need renormalization. However, in the limit 
of vanishing quark masses the (formal continuum) QCD Lagrangian is invariant 
under SU(2)yX SU(2)a flavor symmetry transformations. This leads to non- 
linear relations between the currents called current algebra, from which one 
concludes that no renormalization is necessary (cf. Sect. 6). 

In the regularized theory SU(2)yX SU(2)a is not an exact symmetry but 
is violated by terms of order a. As a consequence there is a finite renormal- 
ization (Meyer and Smith (1983), Martinelli and Yi-Chcng (1983), Groot et al. 
(1984), Gabrielli et al. (1991), Borrclli et al. (1993)) 

(A R )« = Z A A« , 

(V R )l = Z V V« , (6) 

with renormalization constants Z Al Z\? that do not contain any logarithmic (in 
a) or power law divergences and do not depend on any physical scale. Rather 
they are approximated by 

Z A = l + Z A 1) g 2 Q + ... , 

Zv = l + 4 1) .9o + --- , (7) 

3 Of course, decays of hadrons containing b-quarks are more interesting phenomeno- 
logically, but here our emphasis is on the principle of renormalization. 
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for small go- 

On the non-perturbative level these renormalizations can be fixed by current 
algebra relations (Bochicchio et al. (1985), Maiani and Martinelli (1986), Liischer 
et al. (1997 I)) as will be explained in section 6. 

1.3 Scale dependent renormalization 

a) Short distance parameters of QCD. As we take the relevant length 
scales in correlation functions to be small or take the energy scale in scattering 
processes to be high, QCD becomes a theory of weakly coupled quarks and 
gluons. The strength of the interaction may be measured for instance by the 
ratio of the production rate of three jets to the rate for two jets in high energy 
e + e~ collisions, 

a(q) oc ^""""-f , i = (Pe- +Pe + f » lOGcV 2 . (8) 
cr(e+ e — > q q) 

We observe the following points. 

— The perturbative renormalization group tells us that a(q) decreases loga- 
rithmically with growing energy q. In other words the renormalization from 
the bare coupling to a renormalized one is logarithmically scale dependent. 

— Different definitions of a are possible; but with increasing energy, a depends 
less and less on the definition (or the process). 

— In the same way, running quark masses fa acquire a precise meaning at high 
energies. 

— Using a suitable definition (scheme), the g-depcndence of a and fn can be 
determined non-perturbatively and at high energies the short distance pa- 
rameters a and fn can be converted to any other scheme using perturbation 
theory in a. 

Explaining these points in detail is the main objective of my lectures. For now 
we proceed to give a second example of scale dependent renormalization. 

b) Weak hadronic matrix elements of 4-quark operators. Another exam- 
ple of scale dependent renormalization is the 4-fermion operator, As=2 , which 
changes strangeness by two units. It originates from weak interactions after in- 
tegrating out the fields that have high masses. It describes the famous mixing 
in the neutral Kaon system through the matrix element 

(K°|0 As=2 (M)|K°) . 
Here the operator renormalized at energy scale \i is given by 

0^ =2 ( M ) - Z^= 2 ( M a, 5o ) < V^7>d V^7>d + £ z$ s r^ A ~$ s r^ d 

I j=S,P,V,A,T 
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r s = 1, r P 

*j = (fo)> 



75) , 

75, ■■■,r T = cr l 



Zy = —ZA 



(9) 



where I have indicated the flavor index of the quarks explicitly. A mixing of the 
leading bare operator, fp s j^dfp s J^i'd, with operators of different chirality is 
again possible since the lattice theory does not have an exact chiral symmetry 
for finite values of the lattice spacing. The mixing coefficients zj may be fixed 
non-perturbatively by current algebra (Aoki et al. (1997)). Afterwards, the over- 
all scale dependent renormalization has to be treated in the same way as the 
renormalization of the coupling. 

1.4 Irrelevant operators 

A last category of renormalization is associated with the removal of lattice dis- 
cretization errors such as the 0(a)-term in (4). Following Symanzik's improve- 
ment program, this can be achieved order by order in the lattice spacing by 
adding irrelevant operators, i.e. operators of dimension larger than four, to the 
lattice Lagrangian (Symanzik (1982-83)). The coefficients of these operators are 
easily determined at tree level of perturbation theory, but in general they need 
to be renormalized. 

In this subject significant progress has been made recently as reviewed by 
Lepage (1996), Sommer (1997). In particular the latter reference is concerned 
with non-perturbative Symanzik improvement and uses a notation consistent 
with the one of these lectures. It will become evident in later sections that 
improvement is very important for the progress in lattice QCD. 

Note also the alternative approach of removing lattice artifacts order by 
order in the coupling constant but non-perturbatively in the lattice spacing a as 
recently reviewed by Niedermayer (1997). 

2 The problem of scale dependent renormalization 

Let us investigate the extraction of short distance parameters (Section 1.3a) in 
more detail. First we analyze the conventional way of obtaining a from exper- 
iments. Then we explain how one can compute a at large energy scales using 
lattice QCD. 

2.1 The extraction of a from experiments 

One considers experimental observables Oi depending on an overall energy scale q 
and possibly some additional kinematical variables denoted by y. The observables 
can be computed in a perturbative series which is usually written in terms of 
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the MS coupling a-^g, 4 

Oi(q,y) = a m (q)+A i (y)a^(q) + ... . (10) 

For example Oi may be constructed from jet cross sections and y may be related 
to the details of the definition of a jet. 

The renormalization group describes the energy dependence of a in a general 
scheme (a = g 2 /(47r)), 

where the /3-function has an asymptotic expansion 

m s ~ -f{bo + g 2 b 1 + ...} , 

^0 = 74^(11-1^) . 

bi= (^(102-f^) , (12) 

with higher order coefficients &j, i > 1 that depend on the scheme. (12) entails 
the aforementioned property of asymptotic freedom: at energies that are high 
enough for (12) to be applicable and for a number of quark flavors, iVf, that 
is not too large, a decreases with increasing energy as indicated in Fig. 1. The 
asymptotic solution of (11) is given by 

_ 2 9 ^co 1 _ &iln[ln( g 2 M 2 )] / {ln[ln(g 2 M 2 )]} 2 \ 

9 ~ 6 hi(<Z 2 M 2 ) 63[ln(<7 2 /yl 2 )] 2 + \ [\n{q> / A^f ^ 

with A an integration constant which is different in each scheme. 

We note that - neglecting experimental uncertainties - ct^jg extracted in this 
way is obtained with a precision given by the terms that are left out in (10). 
In addition to a 3 -terms, there are non-perturbative contributions which may 
originate from "renormalons" , "condensates" (the two possibly being related) , 
"instantons" or - most importantly - may have an origin that no physicist has 
yet uncovered. Empirically, one observes that values of a-^g determined at differ- 
ent energies and evolved to a common reference point using the renormalization 
group equation (11) including &2 agree rather well with each other; the aforemen- 
tioned uncertainties are apparently not very large. Nevertheless, determinations 
of a are limited in precision because of these uncertainties and in particular if 
there was a significant discrepancy between a determined at different energies 
one would not be able to say whether this was due to the terms left out in (10) or 
was due to terms missing in the Standard Model Lagrangian, eg. an additional 
strongly interacting matter field. 

It is an obvious possibility and at the same time a challenge for lattice QCD 
to achieve a determination of a in one (non-perturbatively) well defined scheme 

4 We can always arrange the definition of the observables such that they start with a 
term a. 
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a 



i i i i i i i i i i i i i i i i i i i 



i i i i I i i i i I i i i i I i i i i I i i i 
log(q) 



Fig. 1. Running of a in a definite scheme. 



and evolve this coupling to high energies. There it may be used to compute jet 
cross sections and compare to high energy experiments to test the agreement 
between theory and experiment. Since in the lattice regularization QCD is nat- 
urally renormalized through the hadron spectrum, such a calculation provides 
the connection between low energies and high energies, verifying that one and 
the same theory describes both the hadron spectrum and the properties of jets. 

Note. A dis-satisfying property of a^jg is that it is only defined in a pertur- 
bative framework; strictly speaking there is no meaning of phrases like "non- 
perturbative corrections" in the extraction of a-^g from experiments. The way 
that I have written (10) suggests immediately what should be done instead. 
An observable d itself may be taken as a definition of a - of course with 
due care. Such schemes called physical schemes are defined without ambigui- 
ties. This is what will be done below for observables that are easily handled in 
MC-simulations of QCD. For an additional example see Grunberg (1984). 



2.2 Reaching large scales in lattice QCD 

Let us simplify the discussion and restrict ourselves to the pure Yang-Mills theory 
without matter fields in this section. A natural candidate for a non-perturbative 
definition of a is the following. Consider a quark and an anti-quark separated 
by a distance r and in the limit of infinite mass. They feel a force F(r), the 
derivative of the static potential V(r), which can be computed from Wilson 
loops (see e.g. Montvay and Minister (1994)). A physical coupling is defined as 

«qq(?) = &r*nr) > 1 = V'. C F = 4/3 . (14) 
It is related to the MS coupling by 

a qq = a us + c f Sqqa lrs + • • • ' ( 15 ) 
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where both couplings are taken at the same energy scale and the coefficients in 

their perturbative relation are pure numbers. The 1-loop coefficient, c^ Sq ^, also 
determines the ratio of the yl-parameters vs. 

A m /A m = cM-cf l «y(8nb )) . (16) 
Note that a q q is a renormalized coupling defined in continuum QCD. 

Problem. If we want to achieve what was proposed in the previous subsection, 
the following criteria must be met. 

— Compute a q q(q) at energy scales of q ~ 10 GeV or higher in order to be able 
to make the connection to other schemes with controlled perturbative errors. 

— Keep the energy scale q removed from the cutoff a" 1 to avoid large dis- 
cretization effects and to be able to extrapolate to the continuum limit. 

— Of course, only a finite system can be simulated by MC. To avoid finite 
size effects one must keep the box size L large compared to the confinement 
scale K^ 1 / 2 to avoid finite size effects. Here, K denotes the string tension, 
K = ]hn r _ 00 F(r). 

These conditions are summarized by 

L y> oAn >y \ ~ locbv >>a < (17) 

which means that one must perform a MC-computation of an TV 4 lattice with 
N = L/a >• 25. It is at present impossible to perform such a computation. The 
origin of this problem is simply that the extraction of short distance parameters 
requires that one covers physical scales that are quite disparate. To cover these 
scales in one simulation requires a very fine resolution, which is too demanding 
for a MC-calculation. 

Of course, one may attempt to compromise in various ways. E.g. one may 
perform phcnomenological corrections for lattice artifacts, keep 1/q ~ a and at 
the same time reduce the value of q compared to what I quoted in (17). Calcu- 
lations of a q q along these lines have been performed in the Yang-Mills theory 
(Michael (1992), Booth et al. (1992), Bali and Schilling (1993)). It is difficult to 
estimate the uncertainties due to the approximations that are necessary in this 
approach. 



Solution. Fortunately these compromises can be avoided altogether (Liischer, 
Weisz and Wolff (1991)). The solution to the problem is to identify the two 
physical scales, above, 

q = l/L . (18) 

In other words, one takes a finite size effect as the physical observable. The 
evolution of the coupling with q can then be computed in several steps, changing 
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q by factors of order 2 in each step. In this way, no large scale ratios appear and 
discretization errors are small for L/a^> 1. 

For illustration, we modify the definition of a q q(<z) to fit into this class of 
finite volume couplings. Consider the Yang-Mills theory on a T x L 3 - torus 
with T 3> L. 5 The finite volume coupling, 

a m (q) ee k{r 2 F(r, L)} r=L/i , q = 1/L , (19) 

can again be related to the MS coupling perturbatively, 

a q ^a m + cf E ^a 2 m + ... . (20) 

This relation may come as a surprise since it relates a small volume quantity 
to an infinite volume one. Remember, however, that once the bare coupling and 
masses are eliminated there are no free parameters. Renormalized couplings in 
finite volume and couplings in infinite volume are in one-to-one correspondence. 
When they are small they can be related by perturbation theory. In particular, 
(16) holds with the obvious modification. 

The complete strategy to compute short distance parameters is summarized 
in Fig. 2. One first rcnormalizes QCD replacing the bare parameters by hadronic 



Lmax = O(ifm) : HS — > SF(q = 1/L max ) 

I 

SF(q = 2/L max ) 
I 



I 

SF(q = 2™/L max ) 
PT: I 

PT 

jet — physics < — Aqcd,M 
Fig. 2. The strategy for a non-perturbative computation of short distance parameters. 



observables. This defines the hadronic scheme (HS) as explained in Sect. 1.1. At 

5 It is well known that perturbation theory in small volumes with periodic boundary 
conditions is complicated by the occurrence of zero modes (Gonzales- Arroyo et al. 
(1983), Liischer (1983)). These can be avoided by choosing twisted periodic boundary 
conditions in space ('t Hooft (1979,1981)), Baal (1983), Liischer and Weisz (1985- 
86)). 
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a low energy scale q = 1/L max this scheme can be related to the finite volume 
scheme denoted by SF in the graph. Within this scheme one then computes the 
scale evolution up to a desired energy q = 2 n /L max . As wc will sec it is no 
problem to choose the number of steps n large enough to be sure that one is 
in the perturbative regime. There perturbation theory (PT) is used to evolve 
further to infinite energy and compute the /1-parameter and the renormalization 
group invariant quark masses. Inserted into perturbative expressions these pro- 
vide predictions for jet cross sections or other high energy observables. In the 
graph all arrows correspond to relations in the continuum; the whole strategy is 
designed such that lattice calculations for these relations can be extrapolated to 
the continuum limit. 

For the practical success of the approach, the finite volume coupling (as well 
as the corresponding quark mass) must satisfy a number of criteria. 

— They should have an easy perturbative expansion, such that the /3-function 
(and r-function, which describes the evolution of the running masses) can 
be computed to sufficient order. 

— They should be easy to calculate in MC (small variance!). 

— Discretization errors must be small to allow for safe extrapolations to the 
continuum limit. 

Careful consideration of the above points led to the introduction of renormal- 
ized coupling and quark mass through the Schrodinger functional (SF) of QCD 
(Luscher et al. (1992), Luscher et al. (1993-94), Sint (1994-95), Jansen et al. 
(1996)). We introduce the SF in the following section. In the Yang-Mills the- 
ory, an alternative finite volume coupling was introduced in G. de Divitiis et al. 
(1994) and studied in detail in G. de Divitiis et al. (1995 I), G. de Divitiis et al. 
(1995 II) . 

The criteria (17) apply quite generally to any scale dependent renormaliza- 
tion, e.g. the one described in Sect. 1.3 b. Although the details of the finite size 
technique have not yet been developed for these cases, the same strategy can be 
applied. This will certainly be the subject of future research. So far, the approach 
has been to search for a "window" where q is high enough to apply PT but not 
too close to a -1 (Martinelli et al. (1994)). An essential advantage of the details 
of the approach of Martinelli et al. (1994) as applied to the renormalization 
of composite quark operators is its simplicity: formulating the renormalization 
conditions in a MOM-scheme, one may use results from perturbation theory in 
infinite volume in the perturbative part of the matching. Since, however, high 
energies q can not be reached in this approach, we will not discuss it further and 
refer to Donini et al. (1995), Oelrich ct al. (1997) for an account of the present 
status and further references, instead. In particular, in the latter reference it 
can be seen, how non-trivial it is to have a "window" where both perturbation 
theory can be applied and lattice artifacts are small. 

Note. (17) has been written for the Yang-Mills theory. In full QCD, finite size 
effects will be more important and one should replace \f~K — > m T , resulting in a 
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3 The Schrodinger functional 

We want to introduce a specific finite volume scheme that fulfills all the require- 
ments explained in the previous section. It is defined from the SF of QCD, which 
we introduce below. For simplicity we restrict the discussion to the pure gauge 
theory except for Sect. 3.7 and Sect. 3.8. Apart from the latter subsections, the 
presentation follows closely Liischer et al. (1992); we refer to this work for further 
details as well as proofs of the properties described below. 



time 




space 

(LxLxL box with periodic b.c.) 



Fig. 3. Illustration of the Schrodinger functional. 



3.1 Definition 

Here, we give a formal definition of the SF in the Yang-Mills theory in continuum 
space-time, noting that a rigorous treatment is possible in the lattice regularized 
theory. 

Space-time is taken to be a cylinder illustrated in Fig. 3. We impose Dirichlet 
boundary conditions for the vector potentials 6 in time, 

A k (x) = i^ a l Xo = ° T , (21) 
feV ' \ C' k (x) at i = i 

where C, C are classical gauge potentials and A A denotes the gauge transform 
of A, 

A A (x ) = yl(x)A fc (x)yl(x)- 1 + ^(x^^x)" 1 , A G SU(N) . (22) 

6 We use anti-hermitian vector potentials. E.g. in the gauge group SU(2), we have 
A fl (x) = A°(i)r° / (2i) , in terms of the Pauli-matrices r a . 
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In space, we impose periodic boundary conditions, 

A k (x + Lk) = A k (x), yl(x + Lk) = A(x) . (23) 

The (Euclidean) partition function with these boundary conditions defines the 
SF, 

Z[C',C] = J T>[A] J B[A}e- Sc ^ , (24) 

Sg[A] = ~ [ d^x tr {F I1V F I1V } , 
z 9o J 

F^v = d^Ay — d v A^ + [Ap, A v \ , 
B[A}= [] dA-{x), D[A]=]\dA(x) . 

Here cL4(x) denotes the Haar measure of SV(N). It is easy to show that the SF 
is a gauge invariant functional of the boundary fields, 

Z[C'°',C n ]=Z[C',C\ , (25) 

where also large gauge transformations are permitted. The invariance under the 
latter is an automatic property of the SF defined on a lattice, while in the 
continuum formulation it is enforced by the integral over A in (24). 

3.2 Quantum mechanical interpretation 

The SF is the quantum mechanical transition amplitude from a state \C) to a 
state \C) after a (Euclidean) time L. To explain the meaning of this statement 
of the SF, we introduce the Schrodinger representation. The Hilbert space con- 
sists of wave-functionals \P[A] which are functionals of the spatial components of 
the vector potentials, A%(x). The canonically conjugate field variables are repre- 
sented by functional derivatives, -EJJ(x) = \ sa^(^) ' an ^ a sca l ar product is given 
by 

W) = J D[AWA}**'[A], D[A] = [] cW£(x) . (26) 

x,/c,a 

The Hamilton operator, 

K = i Ld ^{f^ (x) ^ a(x) + 4^ Ffe " (x) ^ (x) } ' (2?) 

commutes with the projector, IP, onto the physical subspace of the Hilbert space 
(i.e. the space of gauge invariant states), where IP acts as 

TPip[A] = J B[A] ip{A A ] . (28) 
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Finally, each classical gauge field defines a state |C) through 

(C\&) = &[C] . (29) 

After these definitions, the quantum mechanical representation of the SF is given 
by 

Z[C',C] = (C'\e- mT 1P\C) 

OO 

= Y^e- E ^ n [C'\P n [C]* . (30) 

In the lattice formulation, (30) can be derived rigorously and is valid with real 
energy eigenvalues E n . 

3.3 Background field 

A complementary aspect of the SF is that it allows a treatment of QCD in a color 
background field in an unambiguous way. Let us assume that we have a solution 
B of the equations of motion, which satisfies also the boundary conditions (21). 
If, in addition, 

S[A] > S[B] (31) 

for all gauge fields A that are not equal to a gauge transform B n of B, then we 
call B the background field (induced by the boundary conditions). Here, Q{x) 
is a gauge transformation defined for all x in the cylinder and its boundary and 
B n is the corresponding generalization of (22). Background fields B 1 satisfying 
these conditions are known; we will describe a particular family of fields, later. 

Due to (31), fields close to B dominate the path integral for weak coupling 
go and the effective action, 

r[B] = -In Z[C\C] , (32) 
has a regular perturbative expansion, 

r[B] = \r [B] + r 1 [B] + glr 2 [B] + ... , (33) 

So 

r [B}^g 2 Q S[B] . 

Above we have used that due to our assumptions, the background field, B, and 
the boundary values C, C are in one-to-one correspondence and have taken B 
as the argument of r. 
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3.4 Perturbative expansion 

For the construction of the SF-scheme as a renormalization scheme, one needs to 
study the renormalization properties of the functional, Z. Liischcr ct al. (1992) 
have performed a one-loop calculation for arbitrary background field. The cal- 
culation is done in dimensional regularization with 3 — 2e space dimensions and 
one time dimension. One expands the field A in terms of the background field 
and a fluctuation field, q, as 

A„(x) = B„(x) + g q„(x) . (34) 

Then one adds a gauge fixing term ("background field gauge") and the corre- 
sponding Fadeev-Popov term. Of course, care must be taken about the proper 
boundary conditions in all these expressions. Integration over the quantum field 
and the ghost fields then gives 

A [B] = \ In dct A l - In dct A , (35) 

where A\ is the fluctuation operator and A the Fadeev-Popov operator. The 
result can be cast in the form 

A[B] = -hr [B} + 0(1) , (36) 

e— >o £ 

with the important result that the only (for e — > 0) singular term is proportional 
to To- 
After renormalization of the coupling, i.e. the replacement of the bare cou- 
P lin g by g^s via 



the effective action is finite, 



r[B] £=Q = | J- - b Q [hV - ^] \ r Q [B] 



5 MS 

-iC'(0|^i) + C'(0|A>) + O(.gy (38) 



C(0\A) = ±C(s\A) 



C(s\A) = TtA- 



Here, C'(0|^\) is a complicated functional of B, which is not known analytically 
but can be evaluated numerically for specific choices of B. 

The important result of this calculation is that (apart from field independent 
terms that have been dropped everywhere) the SF is finite after eliminating g 
in favor of (^jg. The presence of the boundaries does not introduce any extra 
divergences. In the following subsection we argue that this property is correct in 
general, not just in one-loop approximation. 
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3.5 General renormalization properties 

The relevant question here is whether local quantum field theories formulated on 
space-time manifolds with boundaries develop divergences that are not present in 
the absence of boundaries (periodic boundary conditions or infinite space-time). 
In general the answer is "yes, such additional divergences exist" . In particu- 
lar, Symanzik studied the 4 -theory with SF boundary conditions (Symanzik 
(1981)). In a proof valid to all orders of perturbation theory he was able to show 
that the SF is finite after 

— renormalization of the self-coupling, A, and the mass, m, 

- and the addition of the boundary counter-terms 



In other words, in addition to the standard renormalizations, one has to add 
counter-terms formed by local composite fields integrated over the boundaries. 
One expects that in general, all fields with dimension d < 3 have to be taken 
into account. Already Symanzik conjectured that counter-terms with this prop- 
erty arc sufficient to rcnormalizc the SF of any quantum field theory in four 
dimensions. 

Since this conjecture forms the basis for many applications of the SF to the 
study of renormalization, we note a few points concerning its status. 

— As mentioned, a proof to all orders of perturbation theory exists for the </> 4 
theory, only. 

— There is no gauge invariant local field with d < 3 in the Yang-Mills theory. 
Consequently no additional counter-term is necessary in accordance with the 
1-loop result described in the previous subsection. 

— In the Yang-Mills theory it has been checked also by explicit 2-loop calcula- 
tions (Narayanan and Wolff (1995), Bode (1997)). Numerical, non-perturbative, 
MC simulations (Luscher et al. (1993-94), G. de Divitiis et al. (1995 II) ) give 
further support for its validity. 

— It has been shown to be valid in QCD with quarks to 1-loop (Sint (1994-95)). 

— A straight forward application of power counting in momentum space in 
order to prove the conjecture is not possible due to the missing translation 
invariance. 

Although a general proof is missing, there is little doubt that Symanzik's con- 
jecture is valid in general. Concerning QCD, this puts us into the position to 
give an elegant definition of a rcnormalized coupling in finite volume. 

3.6 Renormalized coupling 

For the definition of a running coupling we need a quantity which depends only 
on one scale. We choose LB such that it depends only on one dimensionless vari- 
able rj. In other words, the strength of the field is scaled as 1/L. The background 




(39) 
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field is assumed to fulfill the requirements of Sect. 3.3. Then, following the above 
discussion, the derivative 

r'[B] = ^r[B] , (40) 

is finite when it is expressed in terms of a renormalized coupling like g^g but r' 
is defined non-perturbatively. From (33) we read off immediately that a properly 
normalized coupling is given by 

f{L) = r'[B]/ r'[B] . (4i) 

Since there is only one length scale L, it is evident that g defined in this way 
runs with L. 

A specific choice for the gauge group SU(3) is the abelian background field 
induced by the boundary values (Liischer et al. (1993-94)) 



C' k = i 4>'i , k = 1,2,3, (42) 





with 



4>i=v-h <P'i = -<t>i- 

(h = ~hv, <A2 = -03 + f 5 (43) 

3 



In this case, the derivatives with respect to ij are to be evaluated at r\ — 0. The 
associated background field, 

B o = 0, B k = [x C' k + (L - x )C k ] /L, fc = l,2,3, (44) 

has a field tensor with non-vanishing components 

Gok - d B k - (C' k - C k )/L, k = 1, 2, 3 . (45) 

It is a constant color-electric field. 



3.7 Quarks 

In the end, the real interest is in the renormalization of QCD and we need to 
consider the SF with quarks. It has been discussed in Sint (1994-95). 

Special care has to be taken in formulating the Dirichlet boundary conditions 
for the quark fields; since the Dirac operator is a first order differential operator, 
the Dirac equation has a unique solution when one half of the components of the 
fcrmion fields are specified on the boundaries. Indeed, a detailed investigation 
shows that the boundary condition 

P+VUo=o = P, P-tp\ X0 =L = p' , P± = i(l± 7 o), (46) 
V;P_L =0 = p, ^P+\ X0 =l = P', (47) 
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lead to a quantum mechanical interpretation analogous to (30). The SF 

Z[C',p',p';C,p,p] = J D[yL]D[^]D[^]e- s ^] (48) 

involves an integration over all fields with the specified boundary values. The 
full action may be written as 

S[A,^] = S G ^] + S F [A,^,i,] 

S F = J d 4 x^(x)[^D^ + m]V(aO (49) 

-J d 3 x$(x)P_Tjj(x)] Xo=0 - j d 3 x$(x)P + iP(x)] Xo=L , 

with Sq as given in (24). In (49) we use standard Euclidean 7-matrices. The 
covariant derivative, Z? M , acts as D fJ ,tp(x) = d^tplx) + A fl (x)^p(x). 

Let us now discuss the renormalization of the SF with quarks. In contrast to 
the pure Yang-Mills theory, gauge invariant composite fields of dimension three 
are present in QCD. Taking into account the boundary conditions one finds (Sint 
(1994-95)) that the counter-terms, 

V>-P-V>|x =o and ipP+ip\ Xo = L , (50) 

have to be added to the action with weight 1 — Z\> to obtain a finite renormalizcd 
functional. These counter-terms are equivalent to a multiplicative renormaliza- 
tion of the boundary values, 

-7-1/2 _/ v -l/2_, /c ,n 

PR = Z h ' p, . .. , p R = Z h ' p . (51) 

It follows that - apart from the renormalization of the coupling and the quark 
mass - no additional renormalization of the SF is necessary for vanishing bound- 
ary values p,...,p'. So, after imposing homogeneous boundary conditions for the 
fcrmion fields, a renormalized coupling may be defined as in the previous sub- 
section. 

As an important aside, we point out that the boundary conditions for the 
fcrmions introduce a gap into the spectrum of the Dirac operator (at least for 
weak couplings). One may hence simulate the lattice SF for vanishing phys- 
ical quark masses. It is then convenient to supplement the definition of the 
renormalizcd coupling by the requirement m = 0. In this way, one defines a 
mass-independent renormalization scheme with simple renormalization group 
equations. In particular, the /3-function remains independent of the quark mass. 



Correlation functions are given in terms of the expectation values of any 
product O of fields, 

(O) =\\\ D[A}D[^}D[^}Oe- s ^^\ , (52) 

I ^ J ) p'=p'=p= P =n 
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evaluated for vanishing boundary values p, . . . , p' . Apart from the gauge field and 
the quark and anti-quark fields integrated over, O may involve the "boundary 
fields" (Liischer et al. (1996)) 

CM = Jr-^ CM 



5p(x)' <5p(x)' 

An application of fcrmionic correlation functions including the boundary fields is 
the definition of the renormalizcd quark mass in the SF scheme to be discussed 
next. 



3.8 Renormalized mass 

Just as in the case of the coupling constant, there is a great freedom in defin- 
ing renormalized quark masses. A natural starting point is the PCAC relation 
which expresses the divergence of the axial current (5) in terms of the associated 
pseudo-scalar density, 

r(x)4(x) 75 lT^(x) , (54) 

via 

d^A^x) = 2mP a (x) . (55) 

This operator identity is easily derived at the classical level (cf. Sect. 6). After 
renormalizing the operators, 

(A K ); = z a a; , 

P£ = Z P P a , (56) 
a renormalized current quark mass may be defined by 

fa = ^-m . (57) 
Z P 

Here, m, is to be taken from (55) inserted into an arbitrary correlation function 
and Z"a can be determined unambiguously as mentioned in Sect. 1.2. Note that m 
does not depend on which correlation function is used because the PCAC relation 
is an operator identity. The definition of fa is completed by supplementing (56) 
with a specific normalization condition for the pseudo-scalar density, fa then 
inherits its scheme- and scale-dependence from the corresponding dependence 
of Pr. Such a normalization condition may be imposed through infinite volume 
correlation functions. Since we want to be able to compute the running mass 
for large energy scales, we do, however, need a finite volume definition. This is 
readily given in terms of correlation functions in the SF. 
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Fig. 4. fp (left) and /i (right) in terms of quark propagators. 

To start with, let us define (isovector) pseudo-scalar fields at the boundary 
of the SF, 

O a = J d 3 u J d 3 v C(u) 75 ^r a C(v), 

O la = /d 3 u| d 3 v C'(u) 7 5^ a C'(v) , (58) 

to build up the correlation functions 

fr{x Q ) = -\{P a {x)O a ) , 

fi = (0 ,a O a ) , (59) 

which are illustrated in Fig. 4. 
We then form the ratio 

Z P = const. yfh/ f P (x)\ X0=L/2 , (60) 

such that the renormalization of the boundary quark fields, (51), cancels out. 
The proportionality constant is to be chosen such that Zp — 1 at tree level. To 
define the scheme completely one needs to further specify the boundary values 
C, C and the boundary conditions for the quark fields in space. These details 
are of no importance, here. 

We rather mention some more basic points about this renormalization scheme. 
Just like in the case of the running coupling, the only physical scale that exists 
in our definitions (57), (60) is the linear dimension of the SF, the length scale, 
L. So the mass fh{L) runs with L. We have already emphasized that g is to 
be evaluated at zero quark mass. It is advantageous to do the same for Zp. 
In this way we define a mass-independent renormalization scheme, with simple 
renormalization group equations. 

By construction, the SF scheme is non-perturbative and independent of a 
specific regularization. For a concrete non-perturbative computation, we do, 
however, need to evaluate the expectation values by a MC-simulation of the 
corresponding lattice theory. We proceed to introduce the lattice formulation of 
the SF. 
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3.9 Lattice formulation 

A detailed knowledge of the form of the lattice action is not required for an 
understanding of the following sections. Nevertheless, we give a definition of the 
SF in lattice regularization. This is done both for completeness and because it 
allows us to obtain a first impression about the size of discretization errors. 

We choose a hyper-cubic Euclidean lattice with spacing a. A gauge field U on 
the lattice is an assignment of a matrix U(x,n) £ SU(iV) to every lattice point 
x and direction n = 0, 1, 2, 3. Quark and anti-quark fields, ip(x) and ip(x), reside 
on the lattice sites and carry Dirac, color and flavor indices as in the continuum. 
To be able to write the quark action in an elegant form it is useful to extend the 
fields, initially defined only inside the SF manifold (cf. Fig. 3) to all times xo by 
"padding" with zeros. In the case of the quark field one sets 

ij)(x) = if x < or x > L, 

and 

P-i>{x)\x =o = P+tp(x)\ Xo =L = 0, 
and similarly for the anti-quark field. Gauge field variables that reside outside 
the manifold are set to 1. 

We may then write the fermionic action as a sum over all space-time points 
without restrictions for the time-coordinate, 

S F [U,i>,^\ =a 4 ^Vi(D + m )V, (61) 

X 

and with the standard Wilson-Dirac operator, 

3 

^lEH^ + ^-^W- (62) 

Here, forward and backward covariant derivatives, 

V^(x) = ±[U(x,^(x + aji)-iP(x)}, (63) 

V^(x) = ±[V0r) - U(x - a£, /*)" V(* - «£)] , (64) 

are used and mo is to be understood as a diagonal matrix in flavor space with 
elements m^. 

The gauge field action Sq is a sum over all oriented plaquettes p on the 
lattice, with the weight factors w(p), and the parallel transporters U{p) around 
P, 

S G [(7] = ^5>(p)tr{l-[/(p)} . (65) 
So V 

The weights w(p) are 1 for plaquettes in the interior and 

( \ _ f 5 c s if P is a spatial plaquette at xo = or x = L, . . 

^ ' \c t Hp is time-like and attached to a boundary plane. ^ ' 
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The choice c s = c t = 1 corresponds to the standard Wilson action. However, 
these parameters can be tuned in order to reduce lattice artifacts, as will be 
briefly discussed below. 

With these ingredients, the path integral representation of the Schrodinger 
functional reads (Sint (1994-95)), 

Z = J V[4,]D[4,]D[U]e- s , S = S F + S G , (67) 
T>[U] = l[dU(x, IJ ,) , 

with the Haar measure dll. 



Boundary conditions and the background field. The boundary conditions 
for the lattice gauge fields may be obtained from the continuum boundary values 
by forming the appropriate parallel transporters from x + ak to x at xq — and 
x = L. For the constant abelian boundary fields C and C that we considered 
before, they are simply 

U(x, fc)U =o = exp(aC fc ), U(x, k)\ Xo=L = cxp(aC' k ), (68) 

for k = 1, 2, 3. All other boundary conditions arc as in the continuum. 

For the case of (42), (43), the boundary conditions (68) lead to a unique (up to 
gauge transformations) minimal action configuration V, the lattice background 
field. It can be expressed in terms of B (44), 

V(x,n) =cxp{aS M (x)} . (69) 



Lattice artifacts. Now we want to get a first impression about the dependence 
of the lattice SF on the value of the lattice spacing. In other words we study 
lattice artifacts. At lowest order in the bare coupling we have, just like in the 
continuum, 

r = \r a [v\ + o((. 9o ) ) , r [v] = g 2 s G [v] . (70) 

Furthermore one easily finds the action for small lattice spacings, 
S G M=[1 + (1- C 0^£ Usin 
3 N 

= 72 [l + (l-c t )f 

We observe: at tree-level of perturbation theory, all linear lattice artifacts are 
removed when one sets c t = 1. Beyond tree- level, one has to tune the coefficient 
Ct as a function of the bare coupling. We will show the effect, when this is done to 



2L 2 



{4>' a - <f> a ) 



0(a 4 



(71) 
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first order in g%, below. Note that the existence of linear 0(a) errors in the Yang- 
Mills theory is special to the SF; they originate from dimension four operators 
FokFok an d FkiF^i which are irrelevant terms (i.e. they carry an explicit factor 
of the lattice spacing) when they are integrated over the surfaces. c s , which can 
be tuned to cancel the effects of F k [F kl , does not appear for the electric field 
that we discussed above. 

Once quark fields are present, there are more irrelevant operators that can 
generate 0(a) effects as discussed in detail in Liischer et al. (1996). Here we 
emphasize a different feature of (71): once the 0(a)-terms are canceled, the 
remaining a-effects are tiny. This special feature of the abelian background field 
is most welcome for the numerical computation of the running coupling; it allows 
for reliable extrapolations to the continuum limit. 



Explicit expression for r'. Let us finally explain that r' is an observable 
that can easily be calculated in a MC simulation. From its definition we find 
immediately 



or] 



|| DMD[^W]e- s j = ^ . (72) 



The derivative || evaluates to the (color 8 component of the) electric field at 
the boundary, 

#f(x) = 4rRe tr h\ 8 U(x,k)U(x + ak,0)U(x + a6, fc)" 1 C/(a;, O)" 1 } 

a 1 V > iro=o 

where As = diag(l, —1/2, —1/2). (A similar expression holds for (E®)'(x)). The 
renormalized coupling is therefore given in terms of the expectation value of a 
local operator; no correlation function is involved. This means that it is easy and 
fast in computer time to evaluate it. It further turns out that a good statistical 
precision is reached with a moderate size statistical ensemble. 



4 The computation of a(q) 

We are now in the position to explain the details of Fig. 2 (Liischer, Weisz and 
Wolff (1991), Liischer et al. (1993-94), Capitani et al. (1997)). The problem has 
been solved in the SU(3) Yang-Mills theory. In the present context, this is of 
course equivalent to the quenched approximation of QCD or the limit of zero 
flavors. We will therefore also refer to results in quenched QCD. 

Our central observable is the step scaling function that describes the scale- 
evolution of the coupling, i.e. moving vertically in Fig. 2. The analogous function 
for the running quark mass will be discussed in the following section. 
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Fig. 5. Schematic plot of the running coupling constructed from the step scaling func- 
tion a. 



4.1 The step scaling function 

We start from a given value of the coupling, u — g 2 (L). When we change the 
length scale by a factor s, the coupling has a value g 2 (sL) = u' . The step scaling 
function, a is then defined as 

a(s,u) = u' . (74) 

The interpretation is obvious. a(s,u) is a discrete /3-function. Its knowledge 
allows for the recursive construction of the running coupling at discrete values 
of the length scale, 

u k =f{s- k L) , (75) 

once a starting value u n — g 2 (L) is specified (cf. Fig. 5). a, which is readily 
expressed as an integral of the /3-function, has a perturbative expansion 

a(s,u) = u + 2b \n(s)u 2 + . . . . (76) 

On a lattice with finite spacing, a, the step scaling function will have an 
additional dependence on the resolution a/L. We define 

E(s,u,a/L) = u' , (77) 

with 

g 2 (L) = u, g 2 (sL) = u' , go fixed, L/a fixed . (78) 

The continuum limit a(s,u) — S(s,u,0) is then reached by performing calcula- 
tions for several different resolutions and extrapolation a/L — > 0. In detail, one 
performs the following steps: 
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1. Choose a lattice with L/a points in each direction. 

2. Tune the bare coupling g such that the renormalized coupling g 2 (L) has 
the value u. 

3. At the same value of go, simulate a lattice with twice the linear size; compute 
u' = g 2 (2L). This determines the lattice step scaling function S(2,u,a/L). 

4. Repeat steps 1.-3. with different resolutions L/a and extrapolate a/L — > 0. 

Note that step 2. takes care of the renormalization and 3. determines the evolu- 
tion of the renormalized coupling. 

Sample numerical results are displayed in Fig. 6. The coupling used is exactly 
the one defined in the previous section and the calculation is done in the theory 
without fcrmions. One observes that the dependence on the resolution is very 
weak, in fact it is not observable with the precision of the data in Fig. 6. We 
now investigate in more detail how the continuum limit of S is reached. As a 
first step, we turn to perturbation theory. 
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u=1.5553 




Fig. 6. Typical example for the lattice step scaling function after 1-loop improvement. 
The continuum limit (circle) is reached by linear extrapolation. 



4.2 Lattice spacing effects in perturbation theory 

Symanzik has investigated the cutoff dependence of field theories in perturbation 
theory (Symanzik (1982-83)). Generalizing his discussion to the present case, one 
concludes that the lattice spacing effects have the expansion 

S(2, u, a/L) — cr(2, u) <-//,-% <■ / , r s ■> 

v 1 ' — = 6 1 (a/L)u + S 2 (a/L)u 2 + ... (79) 
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s n (a/L) a/L ^° J2 e MM!)f(!) + d k AHi )] fe (r) 2 + • ■ • • 

fc=0 

We expect that the continuum limit is reached with corrections 0(a/L) also 
beyond perturbation theory. In this context 0(a/L) summarizes terms that con- 
tain at least one power of a/L and may be modified by logarithmic corrections 
as it is the case in (79). To motivate this expectation recall Sect. 1.4, where we 
explained that lattice artifacts correspond to irrelevant operators 7 , which carry 
explicit factors of the lattice spacing. Of course, an additional a-dependence 
comes from their anomalous dimension, but in an asymptotically free theory 
such as QCD, this just corresponds to a logarithmic (in a) modification. 

t — i — i — i — | — i — i — i — i — r~7 — i — i — i — | — i — i — i — i — 
/ c(»*a/L . 
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Fig. 7. Lattice artifacts at 1-loop order. The circles show 8\(a/L) for the SU(3) 
Yang-Mills theory with 1-loop improvement. The dotted line corresponds to the linear 
piece in a, when only tree-level improvement is used, instead. 
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As mentioned in the previous section, the lattice artifacts may be reduced 
to 0((a/L) 2 ) by canceling the leading irrelevant operators. In the case at hand, 
this is achieved by a proper choice of c t (go)- It is interesting to note, that by 
using the perturbative approximation 

c t ( 50 ) = l + c^ 2 (80) 

one does not only eliminate ei jrl for n = 0, 1 but also the logarithmic terms 
generated at higher orders are reduced, 

e„,„ = 0, e n _i, n = . (81) 

7 For a more precise meaning of this terms one must discuss Symanzik's effective 
theory. We refer the reader to Liischer et al. (1996) for such a discussion. 
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For tree-level improvement, c t (go) = 1, the corresponding statement is e nj „ = 0. 
Heuristically, the latter is easy to understand. Tree-level improvement means 
that the propagators and vertices agree with the continuum ones up to correc- 
tions of order 0(a 2 ). Terms proportional to a can then arise only through a 
linear divergence of the Feynman diagrams. Once this happens, one cannot have 
the maximum number of logarithmic divergences any more; consequently e„ jn 
vanishes. 

To demonstrate further that the abelian field introduced in the previous 
section induces small lattice artifacts, we show 8\(a/L) for the one loop improved 
case. The term that is canceled by the proper choice — —0.089 is shown as a 
dashed line. The left over 0((a/L) 2 )-terms are below the 1% level for couplings 
u < 2 and lattice sizes L/a > 6. We now understand better why the a/ In- 
dependence is so small in Fig. 6. 

From the investigation of lattice spacing effects in perturbation theory one 
expects that one may safely extrapolate to the continuum limit by a fit 

E(2, u, a/L) = cr(2, u) + const, x a/L , (82) 

once one has data with a weak dependence on a/L, like the ones in Fig. 6. Such 
an extrapolation is shown in the figure. 

4.3 The continuum limit universality 
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Fig. 8. Universality test in the SU(2) Yang Mills theory. 

Before proceeding with the extraction of the running coupling, we present 
some further examples of numerical investigations of the approach to the contin- 
uum limit - and its very existence (Liischer et al. (1993-94), G. de Divitiis et al. 
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(1995 II) ). The first example is the step scaling function in the SU(2) Yang-Mills 
theory (G. de Divitiis et al. (1995 II) ). Here we can compare the step scaling 
function obtained with two different lattice actions, one using tree-level 0(a) 
improvement and the other one using c t at ldoop order. (Fig. 8). 

Not only does one observe a substantial reduction of the 0(a)-errors through 
perturbative improvement, but the very agreement of the two calculations when 
extrapolated to a = 0, leaves little doubt that the continuum limit of the SF 
exists and is independent of the lattice action. In turn this also supports the 
statement that the SF is renormalized after the renormalization of the coupling 
constant. 

Turning attention back to the gauge group SU(3), we show the calculation 
of ct(2, u) for a whole series of couplings u in Fig. 9. 
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Fig. 9. Continuum extrapolation of <j(2, u) in the SU(3) Yang-Mills theory. 



4.4 The running of the coupling 

We may now use the continuum step scaling function to compute a series of 
couplings (75). We start at the largest value of the coupling that was covered by 
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Fig. 10. The running coupling in SU(3) Yang-Mills theory. Uncertainties are smaller 
than the size of the symbols. 

the calculation: g 2 = 3.48. This defines the largest value of the box size, L max , 



The series of couplings is then obtained for Lk = 2~ fc L max , k = 0, 1, ... 8. It is 
shown in Fig. 10 translated to a(q) = g 2 (L) / (An) , q = 1/L (We will explain 
below, how one arrives at a GeV-scale in this plot). The range of couplings 
shown in the figure is the range covered in the non-perturbative calculation of 
the step scaling function. Thus no approximations are involved. For comparison, 
the perturbative evolution is shown starting at the smallest value of a that was 
reached. To be precise, 2-loop accuracy here means that we truncate the (3- 
function at 2 loops and integrate the resulting renormalization group equation 
exactly. Thanks to the recent work (Bode (1997), Luscher and Weisz (1995)), 
we can also compare to the 3-loop evolution of the coupling. 

It is surprising that the perturbative evolution is so precise down to very low 
energy scales. This property may of course not be generalized to other schemes, 
in particular not to the MS-scheme, where the /3-function is only defined in 
perturbation theory, anyhow. 

4.5 The low energy scale 

In order to have the coupling as a function of the energy scale in physical units, 
we need to know i max in fm, the first horizontal relation in Fig. 2. In QCD, 



5 2 (£max) = 3.48 . 



(83) 
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this should be done by computing, for example, the product m p L max with m p 
the proton mass and then inserting the experimentally determined value of the 
proton mass. 

At present, results like the ones shown in Fig. 10 are available for the Yang- 
Mills theory, only. Therefore, strictly speaking, there is no experimental observ- 
able to take over the role of the proton mass. As a purely theoretical exercise, 
one could replace the proton mass by a glueball mass; here, we choose a length 
scale, r , derived from the force between static quarks, instead (Sommer(1994)). 
This quantity can be computed with better precision. Also one may argue that 
the static force is less influenced by whether one has dynamical quark loops in 
the theory or not. 




Fig. 11. The dimensionless combination r 2 Fir). The different curves show phenomeno- 
logically successful potential models (Eichten et al. (1980), Martin (1980), Quigg and 
Rosner (1977)). The labels on the top of the graph give the approximate values of the 
r.m.s-radii of the bound states. 



On the theoretical level, ro, has a precise definition. One evaluates the force 
F(r) between an external, static, quark-anti-quark pair as a function of the 
distance r. The radius r is then implicitly defined by 

r 2 F{r)\ r=ra = 1.65 . (84) 

On the other hand, to obtain a phenomenological value for r , one needs to 
assume an approximate validity of potential models for the description of the 
spectra of cc and bb mesons. This is illustrated in Fig. 11. In fact, the value 
1.65 on the r.h.s. of (84) has been chosen to have ro = 0.5 fm from the Cornell 
potential. This is a distance which is well within the range where the observed 
bound states determine (approximately) the phenomenological potential. 
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Fig. 12. Continuum extrapolation of Lmax/ro, using data of Wittig (1995-96), Liischer 
ct al. (1993-94). 

In the following we set r = 0.5 fm, emphasizing that this is mainly for the 
purpose of illustration and should be replaced by a direct experimental observ- 
able once one computes the coupling in full QCD. 

To obtain L max /r from lattice QCD, one picks a certain value of L/a, tunes 
the bare coupling g n such that g 2 — 3.48. At the same value of go one then 
computes the force F(r) on a lattice that is large enough such that finite size 
effects are negligible for the calculation of F(r) and determines tq. Repeating 
the calculation for various values of L/a one may extrapolate the lattice results 
to zero lattice spacing (Fig. 12) and can quote the energies q in GeV, as done in 
Fig. 10. 

4.6 Matching at finite energy 

Following the strategy of Fig. 2, one finally computes the /1-parameter in the SF 
scheme. It may be converted to any other scheme through a 1-loop calculation. 
There is no perturbative error in this relation, as the yl-parameter refers to 
infinite energy, where a is arbitrarily small. 

Nevertheless, in order to clearly explain the problem, we first consider chang- 
ing schemes perturbatively at a finite but large value of the energy. Before writ- 
ing down the perturbative relation between ax and ay where X, Y label the 
schemes, we note that in any scheme, there is an ambiguity in the energy scale q 
used as argument for a. For example in the SF-scheme, we have set q = 1/L, but 
a choice q = n/L would have been possible as well. This suggests immediately 
to allow for the freedom to compare the couplings after a relative energy shift. 
So we introduce a scale factor s in the perturbative relation, 



a Y {sq) = ax{q) + c? Y {s)[ax(q)} 2 +cZ Y {s)[ax(q)} 3 + ... ■ (85) 
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A natural and non-trivial question is now, which scale ratio s is optimal. A 
possible criterion is to choose s such that the available terms in the perturbative 
series (85) are as small as possible. Since the number of available terms in the 
scries is usually low, we concentrate here on the possibility to set the first non- 
trivial term to zero. When available, the higher order one(s) may be used to test 
the success of this procedure. 











Scheme X 


cf MS (1) 


MS (1) 


C2 {SO) 


qq 


-0.0821 


-2.24 


-2.19 


SF 


1.256 


2.775 


0.27 


SF SU(2) 


0.943 


1.411 


0.058 


TP SU(2) 


-0.558 







Table 1. Examples for perturbative coefficients in (85) for N f — 0. 



So we fix s by requiring cf Y (s) — , which is satisfied for s — s with 

so = exp{-cf y (l)/(87T& )} = A x / A Y , (86) 

a relative shift given by the ratio of the /1-parameters in the two schemes. Ex- 
amples taken from the literature (Luscher et al. (1992), Luscher et al. (1993- 
94), Narayanan and Wolff (1995), Bode (1997), Luscher and Weisz (1995), Fish- 
ier (1977), Sint and Sommer (1996), Billoire (1980), Peter (1997)) are listed in 
Table 1. In the case of matching the SF-scheme to MS, the use of so does in- 
deed reduce the 2-loop coefficient considerably. However for the qq-scheme so is 
close to one and the 2-loop coefficient remains quite big. Not too surprisingly, 
no universal success of (86) is seen. 

A non-perturbative test of the perturbative matching has been carried out 
by G. de Divitiis et al. (1995 II) in the SU(2) Yang-Mills theory, where the SF- 
scheme was related to a different finite volume scheme, called TP. 8 The matching 
coefficient for this case is also listed in Table 1. Non-perturbatively the matching 
was computed as follows. 

— For fixed L/a, the bare coupling was tuned such that gg F (L) = 2.0778 (or 
equivalently asF(9 = !/■£) = 0.1653). 

— At the same bare coupling g^ P (L) was computed. 

— These steps were repeated for a range of a/L and the results for g^ P (L) were 
extrapolated to the continuum. 

The result is shown in Fig. 13. 

8 For the definition of the TP-scheme we refer the reader to the literature (G. de 
Divitiis et al. (1994), G. de Divitiis et al. (1995 I)). 
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Fig. 13. Non-perturbative test of perturbative matching. 

We observe that a naive application of the 1-loop formula with s = 1 falls far 
short of the non-perturbative number (the point with error bar), while inserting 
s = s gives a perturbative estimate which is close to the true answer. Indeed, 
the left over difference is roughly of a magnitude a 3 . 

Nevertheless, without the non-perturbative result, the error inherent in the 
perturbative matching is rather difficult to estimate. For this reason it is very 
attractive to perform the matching at infinite energies, i.e. through the A- 
parameters, where no perturbative error remains. 



4.7 The A parameter of quenched QCD 



We first note that the 4-paramctcr in a given scheme is just the integration 
constant in the solution of the renormalization group equation. This is expressed 
by the exact relation 



l 

f3(x) 



1 

box 3 



bi 
Wx 



(87) 



We may evaluate this expression for the last few data points in Fig. 10 using 
the 3-loop approximation to the /3-function in the SF-scheme. The resulting A- 
values arc essentially independent of the starting point, since the data follow the 
perturbative running very accurately. This excludes a sizeable contribution to 
the /3-function beyond 3-loops and indeed, a typical estimate of a 4-loop term in 
the /3-function would change the value of A by a tiny amount. The corresponding 
uncertainty can be neglected compared to the statistical errors. 

After converting to the MS-scheme one arrives at the result (Capitani et al. 
(1997)) 



A^l 

MS 



251 ± 21McV 



(88) 



where the label (°) reminds us that this number was obtained with zero quark 
flavors, i.e. in the Yang-Mills theory. Since this is not the physical theory, one 
must also remember that the overall scale of the theory was set by putting 
r = 0.5 fm. We emphasize that the error in (88) sums up all errors including the 
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extrapolations to the continuum limit that were done in the various intermediate 
steps. 

4.8 The use of bare couplings 

As mentioned before, the recursive finite size technique has not yet been applied 
to QCD with quarks. Instead, a^jg nas been estimated through lattice gauge 
theories by using a short cut, namely the relation between the bare coupling of 
the lattice theory and the MS-coupling at a physical momentum scale which is of 
the order of the inverse lattice spacing that corresponds to the bare coupling (El- 
Khadra et al. (1992)). Without going too much into details, we want to discuss 
this approach, its merits and its shortcomings, here. The emphasis is on the 
principle and not on the applications, which can be found in J. Shigemitsu 
(1996). So, although the main point is to be able to include quarks, we set 
Nf — in the discussion; more is known in this case! 

The method simply requires that one computes one dimensionful experimen- 
tal observable in lattice QCD at a certain value of the bare coupling go. A popular 
choice for this is a mass splitting in the T-system (Davies (1997)). Using as input 
the experimental mass splitting one determines the lattice spacing in physical 
units. 

Next one may attempt to use the perturbative relation, 



to get an estimate for cx-^g. Here we have already inserted a scale shift Sq (cf. 
Sect. 4.6). Without this scale shift, the 1-loop and 2-loop coefficients in the above 
equation would be very large. In turn this means that the shift, 



is enormous. Furthermore, the series (89) does not look very healthy even after 
employing s . Such a behavior of power expansions in a has also been observed 
for other quantities (Lepage and Mackenzie (1993)). One concludes that ao is a 
bad expansion parameter for perturbative estimates. 

The origin of this problem appears to be a large renormalization between the 
bare coupling and general observables defined at the scale of the lattice cutoff 
1/a. Assuming this large renormalization to be roughly universal, one can cure 
the problem by inserting the non-perturbative (MC) values of a short distance 
observable (Parisi (1981), Lepage and Mackenzie (1993)), the obvious candidate 
being 



"Mslsoa X ) = a + 4A5a 3 + O(a^) + 0(a) , a = 9o/(^) 



(89) 



s = 28.8 , 



(90) 



P = ±(trU(p)) . 



(91) 



In detail, due to the perturbative expansion, 



ln(P) = a + 3.373^ + 17.70ag 



+ ... 



(92) 
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we may define an improved bare coupling, 

which appears to have a regular perturbative relation to a^jg, 

a m {s a- v ) =a n + 0.614o£ + 0(a£) + 0(a) . (94) 

Of course, the point of the exercise is to insert the average, P, obtained in the 
MC calculation into (93). Afterwards one only needs to use the (seemingly) well 
behaved expansion (94). One can construct many other improved bare couplings 
but the assumption is that the aforementioned large renormalization of the bare 
coupling is roughly universal and the details do not matter too much. 

On the one hand, the advantages of (94) are obvious: i) one only needs 
the calculation of a hadronic scale and ii) the 2-loop relation to a^jg is known 
(for n/ = 0). On the other hand, how was the problem of scale dependent 
renormalization (Sect. 2.2) solved? It was not! To remind us, the general problem 
is to reach large energy scales, where perturbation theory may be used in a 
controlled way. In the present context this would require to compute with a series 
of lattice spacings for which a^jg(so a_1 ) is both small and changes appreciably. 
The required lattice sizes would then be too large to perform the calculation. 
Therefore one must assume that the error terms in (94) are small. A particular 
worry is that one may not take the continuum limit - due to the very nature 
of (94), which says that a runs with the lattice spacing. This means that it is 
impossible to disentangle the 0(a) and the O(a^) errors. 

We briefly demonstrate now that this last worry is justified in practice. For 
this purpose we consider the SU(2) Yang-Mills theory, where ccsf was computed 
non-perturbatively and in the continuum limit, as a function of the energy scale 
q in units of ro (G. dc Divitiis et al. (1995 II) ). The results of this computation 
are shown as points with error bars in Fig. 14. We may now compare them to 
the estimate in terms of the improved bare coupling, 

a SF (q) = a n + 0.2'Slal , q = Sofl" 1 , s = 1.871 , (95) 

where the only inputs needed are P as well as the value of r /a since qr = 
s r /a. These estimates are given as circles in the figure. 

In general, and in particular for large values of qro, the agreement is rather 
good. However, for the lower values of qro, significant differences are present, 
which are far underestimated by a perturbative error term a 4 . 

What does this teach us about the method as applied in full QCD? To this 
end, we note that the lattice spacings that are used in the applications of im- 
proved bare couplings in full QCD calculations, correspond to q * r < 15. This 
is the range where we saw significant deviations in our test. In light of this it 
appears to us that the errors that are usually quoted for a^g using this method 
are underestimated. It is encouraging, though, that the values which are ob- 
tained in this way compare well with those extracted from experiments using 
other methods (J. Shigemitsu (1996)). 



(93) 
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Fig. 14. Test of an improved bare coupling in the SU(2) Yang-Mills theory. 



5 Renormalization group invariant quark mass 

The computation of running quark masses and the renormalization group invari- 
ant (RGI) quark mass (Capitani et al. (1997)) proceeds in complete analog to 
the computation of a(q). Since we are using a mass-independent renormaliza- 
tion scheme (cf. Sect. 3.8), the renormalization (and thus the scale dependence) 
is independent of the flavor of the quark. When we consider "the" running mass 
below, any one flavor can be envisaged; the scale dependence is the same for all 
of them. 

The renormalization group equation for the coupling (11) is now accompanied 
by one describing the scale dependence of the mass, 

9^ = r(g) , (96) 

where r has an asymptotic expansion 

r(g) ~f {d + fd, + ...}, d = 8/(4tt) 2 , (97) 

with higher order coefficients di, i > which depend on the scheme. 

Similarly to the /1-parameter, we may define a renormalization group invari- 
ant quark mass, M, by the asymptotic behavior of m, 

M = lim m{2bog 2 r do/2b *> ■ (98) 

It is easy to show that M does not depend on the renormalization scheme. It 
can be computed in the SF-scheme and used afterwards to obtain the running 
mass in any other scheme by inserting the proper [}- and r-functions in the 
renormalization group equations. 
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Fig. 15. The step scaling function for the quark mass. 

To compute the scale evolution of the mass non-perturbatively, we introduce 
a new step scaling function, 

cr P = Zp(2L)/Z P {L) . (99) 

The definition of the corresponding lattice step scaling function and the extrap- 
olation to the continuum is completely analogous to the case of a. The only 
additional point to note is that one needs to keep the quark mass zero through- 
out the calculation. This is achieved by tuning the bare mass in the lattice action 
such that the PCAC mass (55) vanishes. At least in the quenched approxima- 
tion, which has been used so far, this turns out to be rather easy (Liischer et al. 
(1997 II)). 

First results for erp (extrapolated to the continuum) have been obtained 
recently (Capitani et al. (1997)). They are displayed in Fig. 15. 
Applying <jp and a recursively one then obtains the series, 

m(2- k L max )/fn{2L max ) , fc = 0,l,... , (100) 

up to a largest value of k, which corresponds to the smallest g that was con- 
sidered in Fig. 15. From there on, the perturbative 2-loop approximation to the 
r-function and 3-loop approximation to the /3-function (in the SF-scheme) may 
be used to integrate the renormalization group equations to infinite energy, or 
equivalently to g = 0. The result is the renormalization group invariant mass, 

M = m(2M 2 )- d "/ 2b «cxp|-^d 5 [lg-^]| . (101) 

In this way, one is finally able to express the running mass fn in units of the 
renormalization group invariant mass, M, as shown in Fig. 16. M has the same 
value in all renormalization schemes, in contrast to the running mass fn. 
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Fig. 16. The running quark mass as a function of q = 1/L. 

The perturbative evolution is again very accurate down to low energy scales. 
Of course, this result may not be generalized to running masses in other schemes. 
Rather the running has to be investigated in each scheme separately. 

The point at lowest energy in Fig. 16 corresponds to 

M/m= 1.18(2) at L = 2i max . (102) 

Remembering the very definition of the renormalized mass (57), one can use 
this result to relate the renormalization group invariant mass mass and the bare 
current quark mass m on the lattice through 

M = mx 1.18(2) x Z A (g )/Zp(g , 2L max /a) . (103) 

In this last step, one should insert the bare current quark mass, e.g. of the strange 
quark, and extrapolate the result to the continuum limit. This analysis has not 
been finished yet but results including this last step are to be expected, soon. 
To date, the one-loop approximation for the renormalization of the quark mass 
(i.e. an approach similar to what was discussed for the coupling in Sect. 4.8) has 
been used to obtain numbers for the strange quark mass in the MS-scheme. The 
status of these determinations was recently reviewed by T. Bhattacharya and R. 
Gupta (1997). 

6 Chiral symmetry, normalization of currents and 
O ( a) -improvement 

In this section we discuss two renormalization problems that are of quite dif- 
ferent nature. The first one is the renormalization of irrelevant operators, that 
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are of interest in the systematic 0(a) improvement of Wilson's lattice QCD as 
mentioned in Sect. 1.4. The second one is the finite normalization of isovector 
currents (cf. Sect. 1.2). They are discussed together, here, because - at least to a 
large extent - they can be treated with a proper application of chiral Ward iden- 
tities. The possibility to use chiral Ward identities to normalize the currents has 
first been discussed by Bochicchio et al. (1985), Maiani and Martinelli (1986). 
Earlier numerical applications can be found in Martinelli et al. (1993), Pacicllo 
et al. (1994), Henty et al. (1995) and a complete calculation is described be- 
low (Liischer et al. (1997 I)). We also sketch the application of chiral Ward iden- 
tities in the computation of the 0(a)-improved action and currents (Liischer et 
al. (1996), Liischer and Weisz (1996), Liischer et al. (1997 I)). 

Before going into the details, we would like to convey the rough idea of the 
application of chiral Ward identities. For simplicity we again assume an isospin 
doublet of mass-degenerate quarks. Imagine that we have a regularization of 
QCD which preserves the full SU(2)vxSU(2)a flavor symmetry as it is present 
in the continuum Lagrangian of mass-less QCD. In this theory we can derive 
chiral Ward identities, e.g. in the Euclidean formulation of the theory. These 
then provide exact relations between different correlation functions. Immediate 
consequences of these relations are that the currents (5) do not get renormalizcd 
(Za = Zy = 1) and the quark mass does not have an additive renormalization. 

Lattice QCD does, however, not have the full SU(2)vxSU(2)a flavor symme- 
try for finite values of the lattice spacing and in fact no regularization is known 
that does. Therefore, the Ward identities are not satisfied exactly. We do, how- 
ever, expect that the renormalized correlation functions obey the same Ward 
identities as before - up to 0(a) corrections that vanish in the continuum limit. 
Therefore we may impose those Ward identities for the renormalized currents, 
to fix their normalizations. 

Furthermore, following Symanzik, it suffices to a add a few local irrelevant 
terms to the action and to the currents in order to obtain an improved lattice 
theory, where the continuum limit is approached with corrections of order a 2 . 
The coefficients of these terms can be determined by imposing improvement 
conditions. For example one may require certain chiral Ward identities to be 
valid at finite lattice spacing a. 

6.1 Chiral Ward identities 

For the moment we do not pay attention to a regularization of the theory and 
derive the Ward identities in a formal way. As mentioned above these identities 
would be exact in a regularization that preserves chiral symmetry. To derive the 
Ward identities, one starts from the path integral representation of a correlation 
function and performs the change of integration variables 

= 1>(x) + ie a A (x)Sl^(x) + iey(x)Syi>(x) , 
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V>(z) -» ^(x)e 4 Vl £ A^- e v(z)] 

= ^(jc) + *e£(a;)<^(aO + ie^(x)S^(x) , (104) 
where we have taken e\(x), ey(x) infinitesimal and introduced the variations 

SW>(x) = ±r>(x), 6^(x) = -?(x)±r a , 

S a A ij(x) = ir a 75 V(x), 5 a A ^{x) =^{x) ln \r a . (105) 

The Ward identities then follow from the invariance of the path integral rep- 
resentation of correlation functions with respect to such changes of integration 
variables. They obtain contributions from the variation of the action and the 
variations of the fields in the correlation functions. In Sect. 6.3 we will need the 
variations of the currents, 

S^Vfe) = -it ahc V^x\ 5%V*{x) - -ie^A^x), 

6%Al(x) = -te abc A^x), S a A A^(x) = -ie abc V*(x) . (106) 

They form a closed algebra under these variations. 

Since this is convenient for our applications, we write the Ward identities in 
an integrated form. Let R be a space-time region with smooth boundary dR. 
Suppose Oi n t and O ex t are polynomials in the basic fields localized in the interior 
and exterior of R respectively. The general vector current Ward identity then 
reads 

/ da^x) (V^x)O int O cxt ) = - ((<^0i„t) Oext) , (107) 

JdR 

while for the axial current one obtains 

/ da^x) (A a (x)O int O ext ) = - ((S a A O int ) Oext) (108) 

JdR 

+2m f d 4 x (P a (x)O int O cxt ) . 

JR 

The integration measure da^ (x) points along the outward normal to the surface 
dR and the pseudo-scalar density P a (x) is defined by 

P a (x)^{x) lb \ T a ^{x) . (109) 

We may also write down the precise meaning of the PCAC-relation (55). It 
is (108) in a differential form, 

([VJW - 2mP a (x)] Oext) = , (110) 

where now O ext may have support everywhere but at the point x. 

Going through the same derivation in the lattice regularization, one finds 
equations of essentially the same form as the ones given above, but with addi- 
tional terms (Bochicchio et al. (1985)). At the classical level these terms are of 
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order a. More precisely, in (110) the important additional term originates from 
the variation of the Wilson term, a V'V^V ^ip , and is a local field of dimension 5. 
Such 0(a)-corrections are present in any observable computed on the lattice and 
are no reason for concern. However, as is well known in field theory, such opera- 
tors mix with the ones of lower and equal dimensions when one goes beyond the 
classical approximation. In the present case, the dimension five operator mixes 
amongst others also with d^A^(x) and P a (x). This means that part of the classi- 
cal 0(a)-terms turn into 0(<7q) in the quantum theory. The essential observation 
is now that this mixing can simply be written in the form of a renormalization 
of the terms that are already present in the Ward identities, since all dimension 
three and four operators with the right quantum number are already there. 

We conclude that the identities, which we derived above in a formal manner, 
are valid in the lattice regularization after 

— replacing the bare fields A, V, P and quark mass m n by renormalized ones, 
where one must allow for the most general renormalizations, 

(A R )» = Z A A; , (Vk)° - Z V V£ , 

(PrT = Z P P a , m R = Z m m q , m q = m Q - m c , 

— allowing for the usual 0(a) lattice artifacts. 

Note that the additive quark mass renormalization m c diverges like 0(g§/a) for 
dimensional reasons. 

As a result of this discussion, the formal Ward identities may be used to 
determine the normalizations of the currents. We discuss this in more detail in 
Sect. 6.3 and first explain the general idea how one can use the Ward identities 
to determine improvement coefficients. 

6.2 0(a)-improvement 

We refer the reader to Liischer et al. (1996) for a thorough discussion of 0(a)- 
improvement and to Sommcr (1997) for a review. Here, we only sketch how 
chiral Ward identities may be used to determine improvement coefficients non- 
pert urbatively. 

The form of the improved action and the improved composite fields is deter- 
mined by the symmetries of the lattice action and in addition the equations of 
motion may be used to reduce the set of operators that have to be considered 
(Liischer and Weisz (1985)). For 0(a)-improvement, the improved action con- 
tains only one additional term, which is conveniently chosen as (Shcikholcslami 
and Wohlert (1985)) 

6S = a 5 ^c sw ^(x)ta^F^(x)V(x), (111) 

X 
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with F^ u a lattice approximation to the gluon field strength tensor F^ v and one 
improvement coefficient c sw . The improved and renormalizcd currents may be 
written in the general form 



(V R )» = Z v (l + b v am q ){V* + ac v ±(d„ + #)T£,}, 

T^(x) = i$a^T a iP(x) 
(A R )£ = Z A (l + b A am q ){A; + ac A ±(^ + d;)P a }, 
(P R ) a = Z P (l + b P am q )P a . (112) 



(<9p and d* are the forward and backward lattice derivatives, respectively.) 

Improvement coefficients like c sw and ca are functions of the bare coupling, 
go, and need to be fixed by suitable improvement conditions. One considers pure 
lattice artifacts, i.e. combinations of observables that are known to vanish in 
the continuum limit of the theory. Improvement conditions require these lattice 
artifacts to vanish, thus defining the values of the improvement coefficients as a 
function of the lattice spacing (or equivalcntly as a function of go)- 

In perturbation theory, lattice artifacts can be obtained from any (rcnormal- 
ized) quantity by subtracting its value in the continuum limit. The improvement 
coefficients are unique. 

Beyond perturbation theory, one wants to determine the improvement coef- 
ficients by MC calculations and it requires significant effort to take the contin- 
uum limit. It is therefore advantageous to use lattice artifacts that derive from a 
symmetry of the continuum field theory that is not respected by the lattice reg- 
ularization. One may require rotational invariance of the static potential V(r), 
e.g. 

V(r = (2,2,l)r/3)-^(r=(r,0,0)) = 0, 
or Lorenz invariance, 

[£(p)] 2 -[£(o)] 2 -p 2 = o, 

for the momentum dependence of a one-particle energy E. 

For 0(a)-improvement of QCD it is advantageous to require instead that 
particular chiral Ward identities are valid exactly? 

In somewhat more detail, the determination of c sw and ca is done as follows. 
We define a bare current quark mass, m, viz. 

W r = K + ac^ + 0' r) P'. ,113) 



9 As a consequence of the freedom to choose improvement conditions, the resulting 
values of improvement coefficients such as c sw , ca depend on the exact choices made. 
The corresponding variation of c sw , ca is of order a. There is nothing wrong with 
this unavoidable fact, since an 0(a) variation in the improvement coefficients only 
changes the effects of order a 2 in physical observables computed after improvement. 
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When all improvement coefficients have their proper values, the renormalized 
quark mass, defined by the renormalized PCAC-relation, is related to m by 

Z A (1 + b A am q ) 2 

m R — y 77— r ^--m + 0(a z . (114) 

Zp(l + bpam q ) 

We now choose 3 different versions of (113) by different choices for O cxt and/or 
boundary conditions and obtain 3 different values of m, denoted by mi, m 2 , m 3 . 
Since the prefactor in front of m in (114) is just a numerical factor, we may 
conclude that all m, have to be equal in the improved theory up to errors of 
order a 2 . c sw and ca may therefore be computed by requiring 

mi = m 2 = m 3 . (H5) 

This simple idea has been used to compute c sw and ca as a function of go in 
the quenched approximation (Luscher et al. (1997 II)). In the theory with two 
flavors of dynamical quarks, c sw has been computed in this way (Jansen and 
Sommcr (1997)). The improvement coefficient for the vector current, cy, may 
be computed through a different chiral Ward identity (Guagnelli and Sommer 
(1997)). 



6.3 Normalization of isovector currents 

Although the numerical results, which we will show below, have been obtained 
after 0(a) improvement, the normalization of the currents as it is described, 
here, is applicable in general. Without improvement one just has to remember 
that the error terms are of order a, instead of a 2 . For the following, we set the 
quark mass (as calculated from the PCAC-relation) to zero. 

Normalization condition for the vector current. Since the isospin sym- 
metry of the continuum theory is preserved on the lattice exactly, there exists 
also an exactly conserved vector current. This means that certain specific Ward- 
identities for this current are satisfied exactly and fix it's normalization automat- 
ically. It is, however, convenient to use the improved vector current introduced 
above, which is only conserved up to cutoff effects of order a 2 . Its normalization 
is hence not naturally given and we must impose a normalization condition to 
fix Zy. Our aim in the following is to derive such a condition by studying the 
action of the renormalized isospin charge on states with definite isospin quantum 
numbers. 

The matrix elements that we shall consider are constructed in the SF using 
(the lattice version of) the boundary field products introduced in (58) to create 
initial and final states that transform according to the vector representation of 
the exact isospin symmetry. The correlation function 

/v(*o) = ^^ abC (°' a ^ b o( x )° C ) ( 116 ) 
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can then be interpreted as a matrix element of the renormalized isospin charge 
between such states. The properly normalized charge generates an infinitesimal 
isospin rotation and after some algebra one finds that the correlation function 
must be equal to 

h = -7^(0' a O a ) (117) 

up to corrections of order a 2 . The 0(a) counter-term appearing in the definition 
(112) of the improved vector current does not contribute to the correlation func- 
tion fy(x a ). So if we introduce the analogous correlation function for the bare 
current, 

M*o) = ^]>> ahc (0'<% 6 (x)0 c >, (118) 

X 

it follows from (117) that 

ZvMx ) = fi + 0(a 2 ). (119) 

By evaluating the correlation functions /i and fv(xo) through numerical simu- 
lation one is thus able to compute the normalization factor Zy. 



Normalization condition for the axial current. To derive a normalization 
condition for Z\, we consider (108) (for m = 0) and choose 0[ nt to be the axial 
current at some point y in the interior of R. The resulting identity, 

/ da^x) {A«(x)Al(y)O cxt )=ie abc (V„ c (y)0 C xt} , (120) 

JdR 

is valid for any type of boundary conditions and space-time geometry, but we 
now assume Schrodinger functional boundary conditions as before. A convenient 
choice of the region R is the space-time volume between the hyper-planes at 
x = yo ± t. On the lattice we then obtain the relation 10 

a 3 £ e abc ( [(A R ) a Q ( yo + 1, x) - (A R r Q (yo - t, x)] (A n ) b (y)O eKt ) 

X 

= 2i ({V R ) c {y)O cxt ) + 0(a 2 ) . (121) 

After summing over the spatial components of y, and using the fact that the 
axial charge is conserved at zero quark mass (up to corrections of order a 2 ), 
(121) becomes 

a 6 £ e abc ((A R r Q (x)(A R ) b (y)O cxt ) = a 3 £ i ((V R ) c Q (y)O ext ) + 0(a 2 ), (122) 
x,y y 

10 In the context of 0(a)-improvement it has been important here that the fields in the 
correlation functions are localized at non-zero distances from each other. Since the 
theory is only on-shell improved, one would otherwise not be able to say that the 
error term is of order a 2 (cf. sect. 2 of Liischer et al. (1996)). 
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where xo = j/o + t. Wc now choose the field product O cxt so that the func- 
tion fy(yo) introduced previously appears on the right-hand side of (122). The 
normalization condition for the vector current (119) then allows us to replace 
the correlation function /y (j/o) by f\. In this way a condition for Za is ob- 
tained (Liischer et al. (1997 I)). 



: — i — i — i — | — i — i — i — | — i — i — i — | — i — i — i — | — i — i — i — | — i — : 
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g§ 

Fig. 17. Current normalization constants as a function of the bare coupling (Liischer 
et al. (1997 I)). The dotted line is 1-loop perturbation theory (Gabrielli et al. (1991)) 
and the crosses correspond to a version of 1-loop tadpole improved perturbation the- 
ory (Lepage and Mackenzie (1993)). The full line is a fit to the non-perturbative results. 



Lattice artifacts and results. It is now straightforward to compute Zy,Zp^ 
by MC evaluation of the correlation functions that enter in (119), (122). Before 
showing the results, we emphasize one point that needs to be considered carefully. 
The normalization conditions fix Zy and Z& only up to cutoff effects of order a 2 . 
Depending on the choice of the lattice size, the boundary values of the gauge field 
and the other kinematical parameters that one has, slightly different results for 
Zy and Za are hence obtained. One may try to assign a systematic error to the 
normalization constants by studying these variations in detail, but since there is 
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no general rule as to which choices of the kinematical parameters are considered 
to be reasonable, such error estimates are bound to be rather subjective. 

It is therefore better to deal with this problem by defining the normalization 
constants through a particular normalization condition. The physical matrix 
elements of the renormalized currents that one is interested in must then be 
calculated for a range of lattice spacings so as to be able to extrapolate the data 
to the continuum limit. The results obtained in this way are guaranteed to be 
independent of the chosen normalization condition, because any differences in 
the normalization constants of order a 2 extrapolate to zero together with the 
cutoff effects associated with the matrix elements themselves. 

Note that a "particular normalization condition" means that apart from 
choosing the boundary values and geometry of the SF, one has to keep the 
size of the SF-geometry fixed in physical units, for example 

L/r Q = const. (123) 

As shown in Fig. 17, the current normalizations can be obtained with good 
precision (in the quenched approximation) after taking all of these points into 
account (Liischer et al. (1997 I)). 

Coming back to our motivation Sect. 1.2, the results in Fig. 17 now allow for 
the calculation of matrix elements of the weak currents involving light quarks 
without any perturbative uncertainties. 

7 Summary, Conclusions 

We have shown how QCD needs to be renormalized non-perturbatively in order 
to obtain unambiguous predictions that can be compared with experiments. 
Once a non-perturbative definition and calculational technique is available, it 
is in principle quite simple to perform renormalization non-perturbatively. In 
practice, the problem has to be treated with care. 

The only presently available definition is lattice QCD with MC simulations as 
the calculational tool to get predictions. In this case, straightforward solutions to 
the renormalization problem face a serious difficulty: the theory must be treated 
at various different energy scales simultaneously, which is an extremely hard 
(impossible?) task for a MC simulation. To circumvent this difficulty, Liischer, 
Weisz and Wolff have introduced the recursive finite size technique, where one 
connects low and high energies recursively in small steps. We have shown, how 
this idea can be put into practice in QCD using the Schrodinger functional as 
a second technical tool. In the theory without dynamical quarks, these meth- 
ods have been shown to allow for the computation of short distance parameters 
like yl^jg and the renormalization group invariant quark mass with completely 
controlled errors! From the practical point of view, the non-perturbative renor- 
malization of other quantities, such as the As = 2 operator, still have to be 
investigated, but no new difficulties are expected to appear. It is therefore plau- 
sible that the renormalization problem can be solved for many specific cases - 
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and with good accuracy. It should not remain unmentioned, however, that in 
practice each renormalization problem has to be considered separately. Certain 
problems may turn out to be significantly more difficult to solve than the ones 
discussed in the lectures. 

We also sketched, how Symanzik improvement can be implemented in a 
non-perturbative way, reducing the leading lattice artifacts from linear in a to 
quadratic in a (O(a)-improvement). For light quarks, such a project has already 
been done carried out. As a result, significant progress in lattice QCD is expected 
from the use of 0(a)-improved QCD. 
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